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Abstract 

I consider the Hermitean two-matrix model with a logarithmic potential which is 
associated in the one-matrix case with the Penner model. Using loop equations I find 
an explicit solution of the model at large ./V (or in the spherical approximation) and 
demonstrate that it solves the corresponding Riemann-Hilbert problem. I construct 
the potential of the Kazakov-Migdal model on a D-dimensional lattice, which turns 
out to be a sum of two logarithms as well, whose large-A^ solution is given by the 
same formulas. In the "naive" continuum limit this potential recovers in D < 4 
dimensions the standard scalar theory with quartic self-interaction. I exploit the 
solution to calculate explicitly the pair correlator of gauge fields in the Kazakov- 
Migdal model with the logarithmic potential. 
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1 Introduction 



The Penner model [[[[] is the Hermitean one-matrix model defined by the partition function 

Z 1 = Jd<Pe- NtlV ^ (1.1) 

with the logarithmic potential 

F(0) = -a( log (1 + 0) -</>). (1.2) 

This model has been extensively studied in the modern context of matrix models after the 
observation by Distler and Vafa that the partition function (|1.1|) with the potential (|1.2|) 
coincides with that of the compactified one- dimensional string at a self-dual radius. Some 
generalizations of the Penner model [|3] were considered and solved in the large- N limit. 
Of special interest is the fact H that an external field problem for the Penner model is 



equivalent to the Hermitean one-matrix model ( |1 . 1|) with a general potential V(<p) which 
is analytic at <fi = 0. 

In the present paper I consider the Hermitean two-matrix model 

^# 2 e Ntr ( W2 " V(W - V(W ) (1.3) 



with the Penner potential (1^2). One of the motivations is the recent paper || where the 



two-dimensional minimal conformal field theories have been constructed from the two- 
matrix model (|1.3|) with polynomial potentials as was first conjectured by Douglas ||. 
From this point of view it is interesting to find out what happens for the non-polynomial 
potential (|1.2|) . 

Another motivation is that the two-matrix model ( |1.3|) describes the solution of the 
Kazakov-Migdal |7j model which is defined by the partition function 



J KM 



/ n du,(x) n d<p x e^ Nt - w) ) (L4) 

X,fJ, X 



where the integration over the gauge field U^x) is over the Haar measure on SU(N) at 
each link of a D-dimensional lattice with x labeling its sites. The Kazakov-Migdal model 
is an interesting extension of the matrix models ( |1 . 1| ) and ( |1.3| ) to the multi-dimensional 
case and obviously recovers the two- matrix model if the lattice is just two-points since 
the gauge field can be absorbed by a unitary transformation, say of 4>i. 

While the fact that the models ( |1.3| ) and ( |1.4j ) are equivalent at large N is known || , 
an explicit relation between the potentials V in (|1.3|) and V in (l.4|) has been studied 



recently || . The idea of the DMS approach is first to solve the two- matrix model ( |1.3|) in 
the large- N limit for some choice of V which is practically convenient and then, having the 
solution, to calculate the associated potential V of the Kazakov-Migdal model. To realize 
this program, the loop equation for the one-link correlator was derived whose solution 
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completely determines all the correlators of the model, in particular, the pair correlator 
of the gauge fields. Besides the Gaussian case when both potentials are quadratic, some 
explicit relations between V and V have been studied. 

In the present paper I show that for the Penner potential ( |1.2| ) (in fact for a little bit 
more general one given by Eq. ( |4.1| ) below) the associated potential V of the Kazakov- 
Migdal model is given by the sum of two logarithms as well (Eq. Q4-19Q below) so that 
the model ( |1.4|) with this potential is explicitly solved at large N for any D. The solution 
is conveniently described via the one-matrix model whose potential (Eq. ( [4.9|) below) 
coincides with V of the Kazakov-Migdal model for D — 0. I demonstrate that this 
solution satisfies the master field equation derived by Migdal |L0| and recovers by turning 



the parameters of the potential the solution for the quadratic potential found by Gross |TT[ . 
I study some properties of the Kazakov-Migdal model with this logarithmic potential and 
show that in the "naive" continuum limit it reproduces in D < 4 the standard action of 
the scalar matrix field with quartic self-interaction. I calculate explicitly the correlator of 
two gauge fields U^(x) and UUx) at the same link (x, /i) which completely determines all 
other correlators of the model and show that it might undergo phase transitions at the 
values of parameters where the proper one-matrix model does. 



2 A review of the DMS approach 



Let us define for the Kazakov-Migdal model (|1.4|) the loop average 



%(t±t)) (2-1) 



N \\ - 



'X 




which describes all correlators of powers of cj) x at the same site and the one-link correlator 

1 — (2.2) 

Both functions are analytic on the complex plane with some cut (or cuts) at the real axis 
and obey the following asymptotic expansions 



and 



(2.4) 



Notice that G v \ is symmetric in v and A due to invariance of the Haar measure, dU, under 
transformations U — > IP. 
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The correlator G u \ obeys in the large- iV limit the following equation || 

/ n — : — ^-G w \ = E V G V \ + \G V \ — E v , (2.5) 
JCi 2m v — u 

where the contour C\ encircles counterclockwise the cut (or cuts) of the function G v \ (the 
same as the support of p) and 

V'{uj) = V\u) - (2D - l)F(w) . (2.6) 

The function 

oo 

F(u) = £ F n" n (2-T) 

n=0 

is the one which appears in the expansion of the pair correlator of the gauge fields 



where the averaging is only w.r.t. U while <j) and x play the role of external fields. As was 
proposed in Refs. |L0], the following formula holds at iV = oo: 

tr \ 00 tr 

^(^X^)) ff =£^ (2-9) 

where t a (a = 1, . . . , N 2 — 1) stand for the generators of the SU(N) which are normalized 
by 

J2(t a ) ij (t a ) kl = N5 a 5 kj - 5 ij 5 kl . (2.10) 



The structure of the r.h.s. of Eq. ( |2.9|) can be understood analyzing the power series 
expansion of the Itzykson-Zuber correlator in <p and x 

£ = EC f | (*V) , (2.1D 

where the coefficients C nm depend only on the distribution of eigenvalues of and x- This 
formula is due to the invariance of the Haar measure dU, which forbids terms like tr <px, 
as well as to the fact that only the traceless parts of <fi and x contribute to the integral in 
Eq. ( 12.111) . Multiplying Eq. ( p. 11] ) by tit b x and using the completeness condition ( |2.10|) , 
one identifies F n with 

F n = ECn, m (^X m+1 - |x|x m ) , for n>l. (2.12) 
Eq. (|2.6j) arise after the multiplication of Eq. (|2.9|) by tr (t a j(v — (j))) which gives 



- (—UxU^) ) = - (^.) + - (— ) - U - F(</>j) . (2.13) 



u 



Since Fq does not enter this formula, one can always choose 
to provide 



n=l 



tr 

iV 







(2.14) 



(2.15) 



-F{4>) j 

For the even potential V(4>) = V(—(j)) Eq. ( |2.15| ) holds automatically since the distribution 
of eigenvalues is symmetric, p(A) = p(— A), and F(X) = —F(—X). 

With F given by ( 2.14j ), Eq. ( p2 . 1 3[ ) can be rewritten as 

duo F(u) 



where G\(v) is defined by the expansion (|2] 
loop equation (see |TB| for a review) 

r du V\u) 

JCx 2-Kl V — UJ 

which is just a result of the infinitesimal shift 



-E„ 



(2.16) 



Ci 2m' v 

The same quantity appears in the simplest 



E u - 2DG x (v) 



E 



v 



(2.17) 



(2.18) 



of the measure in ( |1.4|) and can be alternatively obtained taking the 1/A term of the 
expansion of Eq. Q2.5|) in A. Inserting ( [2.16 ) into Eq. ( [2.17 ) and using ( [2.6|) , one arrives 
at the following equation for E v 



du V'(uj) 



Ci 2m v 



E„, = E„ 



LO 



;2.i9) 



where 



V'(\) = V'(A) - F(X) . (2.20) 

One can see now that Eq. (|2.15|) can be rewritten using Eqs. fl2.19|) and (|2.20| ) as 

du 

'Ci 

which coincides with the l/(uX) term of the double expansion of Eq. ( |2.5| ) in 1/u and 



(2.21) 



1/A. Therefore, the choice ( ^.14| ) is always compatible with Eq. ( p.5[ ). 

While Eq. ( |2.19| ) coincides with the large- loop equation for the Hermitean one- 
matrix model, the potential V(uj) is, generally speaking, non-polynomial and has sin- 
gularities on the complex plane outside of the cut (or cuts) of E u . Another interesting 
property of the above formulas is that at D = 1/2, which is associated with the Hermitean 
two-matrix model, the last two terms on the r.h.s. of Eq. (|2.6j ) disappear and one gets 
just V(uj) = V{u). For this reason the Schwinger-Dyson equations for the Hermitean 
two-matrix model f\A 



15 1 are equivalent to Eq. 
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3 Exact solution for arbitrary potential 



Let us consider the two-matrix model with an arbitrary potential 

vH = E^ m - (3-i) 

m=l " b 

The solution of loop equations in this case was obtained in Ref. |TJ[ and in a more explicit 



form in Ref. @. The equation Q2.5Q for the potential ( p.l|) reads explicitly 

(\ + E u -V(v))G uX = E v -R x {v) (3.2) 

where R\{y) is given by 

AaM e= - / |^^G. A = E ^ E G„(A)i^-»- 2 (3.3) 
Jc 2 2hiv-uj ^ ^ 

and the contour Ci encircles both singularities of G^x and the pole at w — v. The terms 
on the r.h.s. result from taking the residue at infinity in the contour integral and the 
functions G n (X) are defined by Eq. 

The formal solution to Eq. (|3.2p is 



E v - R x (u) 

G » X ~ X + E u -V'(u) (3 - 4) 

with R\(v) given by ( p.3|) . The functions G n (\) can be expressed via E x using the 
recurrence relation 

G n+ i(X) = I ^G„( W ) - E x G n (X) , G (X) = E x (3.5) 

JCi 2,1X1 A — Id 



which can be obtained expanding Eq. ( |2.5| ) in 1/A. For n = this equation recovers 
Eq. (£TSD. 

To determine E\, I use the equation 

^ 7 V(cj)G uX = XE x -l (3.6) 

which is nothing but the 1/v term of the expansion of Eq. (|2.5| ) in 1/u. Using the 
expansion ( |2.4j) , one can rewrite Eq. (|3.6|) as 

E 9mG m -i{X) = XE x -l. (3.7) 

m>l 



For V(A) being a polynomial of the highest power J, Eq. (|3.7|) contains E\ up to the power 
J and the solution is therefore algebraic. This fact was first noted by Staudacher [[14]] for 
the two-matrix model. An explicit example of a non-polynomial potential which results 
in a quadratic equation for E\ is given below. 

As is proven in Ref. 0: 
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• Equations which appear from the next terms of the 1 / //-expansion of Eq. Q2.5Q are 
automatically satisfied as a consequence of Eqs. ( |3.5| ) and ( |3.7| ). 

• G u \ given by ( |3.4j) is symmetric in v and A for any solution of Eq. (|3.7|) . The sym- 
metry requirement can be used directly to determine E\ alternatively to Eq. (|3.7| ). 

It is worth mentioning the relation to the approach by Migdal JlO| which is based on 
the Riemann-Hilbert method. As is noticed in Ref. [0], Eq. (|2.5|) is equivalent to the 



equation of Ref. [10|. To show this let us define the continuous and discontinuous in v 
parts of G v \ across the cut (cuts) by 

G(v+io)A — G{ 



Disc, G 



Cont, G vX = ' U+W,A ' {U - lV)A (3.8) 

so that for a real A outside of the cut (cuts) Disc,, G v \ coincides with the imaginary part 
and Cont, G,a coincides with the real part. In particular, Disc, .E, = —7rp(v). The 
discontinuous part of Eq. (|2.5|) then reads 



Cont, G u \ = 1 + — ^ (A + \v'{v) - V{u)) Disc, G uX , for v e cut (3.9) 
Tvp{u) V 2 ' 



which coincides with the equation of Ref. WU 



To obtain the solution to Eq. ( |3.9| ) for G v \ versus E v , one notices that for any real u 

1 - G( u+i0 )x _ 1 - Cont, Gux - i Disc, G uX _ A + Re E„ - i Im E„ — V{v) _ ^ 
1 - G(y-ia)\ 1 - Cont, G, A + i Disc, G vX A + Re E v + i Im E v - V'(v) 

since Disc, G v \ cancels at the cut (cuts) due to Eq. ( p.9|) . The analytic function which 
solves this problem is given by (TO, O, [if 



G, A = l-exp{-/ ^^_log(A + ^-V'(a;))i ■ (3-11) 
[ Jd 2m v — u) J 

Taking the residue at u = u, one can rewrite it in the form (|3.4|) with 



^<"> ^— togfl + , E : v( M . (3.12) 



/c*2 2m v — uj \ A — V'{uj 

The l.h.s. is obviously analytic in A, expandable in 1/A, has no discontinuity in v and is 
expandable in powers of v — i.e. has all the properties required for R x (v). 

The formula (|3.11| ) should solve the recurrence relation (|3.5| ). I have checked this 
expanding in 1/A to few lower orders. Analogously, the master field equation 

E X = f ^Llog(X + E u - VH), (3.13) 
JCi 2m 

which is obtained as the 1/v term of Eq. (|3.11|) , should be equivalent to Eq. ( [3.7D . It 
is a matter of practical convenience which form of equations to solve. As was shown 
in Ref. ||, Eqs. (|2.5| ) and fl3.7|) can easily be solved for the quadratic potential V(0). A 
simple solution of Eqs. Q2.5|) , (|3.7|) for the Penner potential is obtained in the next section. 
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4 An explicit solution for Penner potential 



Let us choose the following potential of the two matrix model 

°° ab + c c 

V(0) = -(oft + c) log (6 - 0) - a0 = £ — + 70 (4-1) 

m =2 m0 

where a and 6 are real. This potential recovers the quadratic one in the limit 

a, b — > oo , 5 = A~"^' c ~l (quadratic potential) . (4.2) 



The loop equation fl2.5|) involves 



V» = f^. (4.3) 



The function R\(v) which enters (|3.4|) can then be calculated to give 



j-v / \ /"do; V'(a>) _ ab + c„ 

JC* 2 27TZ Z/ — W — V 

This formula is obtained doing the integration over the contour Ci which encircles both 
singularities of G u \ and the pole at uj = v. For V'{uj) given by ( |4.3| ) the residue at infinity 
vanishes since G w \ decreases as 1/u while the residue at uj = b results in the expression 
on the r.h.s.. One can alternatively derive this formula expanding V'(uj) in uj and taking 
the residue at infinity. 

The function Gb\ on the r.h.s. of Eq. fl4.4|) can be determined from Eq. (|3.6|) to be 



ab + c 

so that the solution to Eq. (|2.5| ) for G v \ versus E\ is 

E _ {a+X)E x -l 

G «= x" +K -?* - (4 - 6) 

^ b—v 

To determine E\ let us utilize the symmetry of G v \ in v and A. From Eq. ( }4.6| ) one 

gets 

(b - v)E v - (a + fe)E fe + 1 (4 = 5) (a + - 1 
1/6 (&_„)(& + £„) -ai/-c afe + c 1 j 

which yields the quadratic equation for E v 

[a + v){b- v)El - ({av + c)(o + v) + (1 + c - 6i/)(6 - z/))-E„ 

= — i>(a + 6) + 6 — c — (afe + c) ((a + b)E b -l). (4.8) 
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I have verified that ( f4.6|) is indeed symmetric in v and A providing Eq. ( (4.8[) is satisfied. 
The formula for C(u, A), the double discontinuity of G v \ both in v and A across the cut 
(cuts), which is manifestly symmetric in v and A is presented in Sect. ^|. 

Eq. ( |4.8| ) exactly coincides with the large- N loop equation for the Hermitean one- 
matrix model with the logarithmic potential 

V{4>) = -(ab + c) log {b - <p) + (ab + c + 1) log (a + <p) - (a + b)<p (4.9) 

for which 

nA) = ^ + i±^. (4.10) 



b — A a + A 

The constant £7& is a free parameter which should be fixed requiring the analytic structure 
of E\ on the complex A-plane. 

It it easy to demonstrate that the solution Q4.6|) satisfies Eq. fl3.9|) . To show this let us 
use the symmetry of (|4.6| ) in v and A and calculate Disc,, G\ v and Cont^ G\ v which read 



1 



DisCj, G 



b-X 



Xu 



v + E x - 



aX+c ' 



Cont„ G 



\{a+v)V'{v)-l 
'' b-X 



v ^Xu 



u + Ex- 



aX+c 



(4.11) 



b-X " 1 ^ 6-A 

These expressions identically satisfy Eq. (|3l9|) and, therefore, Eq. ( |3.10|) . I conclude that 
( |4.6|) solves the Riemenn-Hilbert problem and should coincide with ( |3.11| ) for the potential 
V given by Eq. ( |4.3|) . Analogously, Eq. (|3.13|) should reduce to Eq. ( |4.8|) . 



It is instructive to discuss the one-cut solution to Eq. ( f4.8|) which reads [17 

du V'(u) y/(A-x)(A-y) 



Ex , , . . 

where the ends of the cut, x and y, are determined by the asymptotic conditions 



du 



Ci 2iri 







(u - x)(u - y) 



duo ujV'(uS) 



Ci 2m 



(u - x){u - y) 



2. 



(4.12) 



(4.13) 



For V' given by (|4.10|) the contour integral can easily be calculated taking the residues at 
u = X,b and —a while the residue at infinity vanishes since falls down as One 
gets 



1 i aX + c 

Ex — 

> 



1 + c-bX 



ab + c 



1 



+ 



6-A a + X 
ab + c + 1 



a + X 



yj(-a-x)(-a-y) 



A yj(b-x)(b-y) 
y/(X-x)(X-y). 



(4.14) 



It is worth noting that all the formulas recover the ones for the Gaussian potential in 
the limit ( |4.2|) . Since in the Gaussian case Eq. (|4.13 ) has the solution 

o a b 
b a ' 



y 



X 



(quadratic potential) , 



(4.15) 
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the one-cut solution ( 4.14 ) is always realized for a and b which are big enough for the 
points b and —a to lie outside of the cut [x,y]. 

The situation when ether b or —a touch the end of the cut is associated with a phase 
transition to a more-than-one-cut solution. The one-matrix model with the potential (|4.9|) 
has at least as rich structure as that with the cubic potential since the latter can be 



obtained from ( |4.9| ) taking the limit 

a, b — > oo, (a — b) ~ &3 (cubic potential) (4.16) 

and rescaling the field — > 0&s which can always be done for the one-matrix model. The 
linear in term is finite providing c < bt . It can be explicitly seen that in the limit (|4.16|) 
the one-cut (|4.14j) solution recovers the one for the cubic potential. 

Since V'(A) is known, the function F(X) can be determined from Eq. (|2.20|) to be 



™ = ^. (4,7) 



Now Eq. ( P^BD determines 



V W = £±£ + (2D - D^fi (4.18) 
— A a + A 

which corresponds to the potential 
1/(0) = -(ab + c)\og(b-<p)-(2D-l)(ab + c+l)\og(a + <p) + [(2D-l)b-a]<p. (4.19) 



This formula recovers at D = 1/2 the potential (O) of the two- matrix model and at 
= the potential ( (4.9[) of the associated one-matrix model. In the Gaussian limit (|4.2|) 
the solution found by Gross is reproduced. Analogously, I verified that the next term 



of the 1/6 expansion recovers a perturbation of the model with the quadratic potential 
by a cubic term with the coupling g 3 = a/b 2 . 



5 The "naive" continuum limit 



The potential ( 4.19| ) of the D-dimensional Kazakov-Migdal model admits the following 



"naive" continuum limit. Let us expand V in which yields 



i.2 



+ \{^ + i 2D-l)^±l)^ + Om. (5.1) 
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Now the idea is to choose the couplings a, b and c to provide the canonical scaling 

= e f- 1 $ ) v(<f>) = e D W(^) (5.2) 



which is prescribed, as usual, by the kinetic term. The expansion will be justified choosing 
a ~ b — > oo so that 0/6^0. The linear in term has the continuum limit providing 

~ + 1 - 2D) c - (2D - 1) < e% +1 a (5.3) 



,6 

which can always be satisfied by a proper choice of c. 
The bare mass 

m 2 = % + (2D - 1)- > 2D (5.4) 
o a 

for any .D and the equality sign holds only for D — 1 at a/b — 1. Since the equality is 
needed to cancel the term 2D<p 2 coming from the kinetic energy, one concludes that the 
"naive" continuum limit is possible only for D = 1. 

At D = 1 one has 

ab 

where 

^-^ = M ~ 1 (5.6) 

eb 

is the continuum mass. Now the cubic in term has the continuum limit providing 

a = b = - 2 (5.7) 

which justifies the expansion in 

(5.8) 

Notice that the quartic in term is also finite for the ^-dependence of b given by Eq. ( |5.7| ) 
as well as Eq. ( |5.3|) is satisfied for any c ~ 1. We arrive, therefore, at the one-dimensional 
continuum quartic action 

S[$) = J di(const.$(t) + ^$ 2 (t) + ^$ 2 (t) + ^M&{t) + |$ 4 (t)) (5.9) 

where $ = d t <&. The coefficient in front of the linear in term can be made arbitrary (or 
vanishing) by choosing c. 

The above procedure of taking the "naive" continuum limit works for D < 4 where 
b given by Eq. (|5.7|) is divergent. The only difference is that for D ^ 1 the bare mass 
term remains infinite. One can cancel it, however, for 2 < D < 4 by the standard 
renormalization procedure. The restriction D < 4 is precisely the one where the scalar 
theory with the quartic interaction is renormalizable. 
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As was shown by Gross |fLl j, any solution of the master field equation (|3.13|) for 
an arbitrary potential reproduces in D — 1 the spectral density given by the fermionic 
solution of Ref. 



For this reason there is no doubt that the solution of the previous 
section does. However, an interesting question is how the one-cut solution ( |4.14| ) can 
1 spectral density? The answer is that the potential (5J}) is simple 
1 the factorized structure 



reproduce this D 
enough to provide at D 



vrp(A) = y/2E F - 2W(X) = </f (A + const. )yj(y - A)(A - x) 



(5.10) 



by turning the coupling in front of the linear in <fi term. 



6 The Itzykson— Zuber correlator 

To calculate the pair correlator of U and TP, one should take the discontinuity of (|4.6|) 
both in v and A across the cut (cuts) f9j. For the solution ( |4.6|) one gets 

C ( u i ^) = o / \ / \ \ Disc A Disc f G v\ = ( a ^ , u, A G cut (6.1) 



with 



V(v, A) = -AV + {b- a)Xv{X + v)+ ab{\ 2 + v 2 ) - \v{a 2 + b 2 + 2c + 1) 

+ (A + v){b - ac + be) - c 2 -b 2 + (ab + c) [(a + b)E b - 1] , (6.2) 



where E b depends on the type of the solution of Eq. (|4.8p (one-cut or more-than-one-cut 
solutions). In the Gaussian limit (|4.2j ) when E b — > 1/6, one recovers the result of Ref. ||. 

While the expression (|6.1|) looks complicated, the denominator seems to be positive. 
To see this let us rewrite T){y, A) as 

V{v, A) = (b - u)(a + u)(X- r+(p))(X - r_(i/)) (6.3) 

and 

r±{y) = \{F{u) + V{y)) ± iitp{v) . (6.4) 

Since p(v) is real for v G cut, the roots are complex except for the values of parameters 
when p{y) vanishes for some v G cut. This might happen for the cases of phase transitions 
in the one-matrix model with the potential ( (4.9| ). The simplest one is from the one-cut 
solution to a two-cut solution. It is an algebraic problem to study whether the roots ( J6.4Q 
always lie outside of the cut (cuts). 

These phase transitions does not change, however, the functional structure of C{y, A) 
given by Eqs. (|6.1|) , Q6.2p — only the constant E b changes. For this reason the Kazakov- 
Migdal model with the logarithmic potential (|4.19 ) seems to be too simple to exhibit a 
non-trivial critical behavior at D > 1. 
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